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We compute the flavorless running coupling constant of QCD from the 
three gluon vertex in the (regularisation independent) momentum subtrac- 
tion renormalisation scheme. This is performed on the lattice with high 
statistics. The expected color dependence of the Green functions is verified. 



On 



\ There are significant 0(a 2 fi 2 ) effects which can be consistently controlled. 

Scaling is demonstrated when the renormalisation scale is varied between 
2.1 GeV and 3.85 GeV. Scaling when the lattice spacing is varied is very 

^ (c) 

well satisfied. The resulting flavorless conventional two loop Ai^L is esti- 
mated to be, respectively for the MOM and MOM scheme, 361(6) MeV 

(3) 

^ . and 345(6) MeV , while the three loop results are, depending on = 

(412-59-^±6MeV) and A^L = (382 -46 ^=±5 MeV). A preliminary 

computation of (32 in the MOM scheme leads to A^ = 303(5) a 1 qjq^ 



MeV. 
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The non-perturbative calculation of the running coupling constant of QCD is 
certainly one very important problem. This program has been performed using 
the Schrodinger functional |l|], the heavy quark potential the Wilson loop 

[|J, the Polyakov loop || and the three gluon coupling ||. The latter method is 
the one we will follow in the present letter. The principle of the method is quite 
simple since it consists in following the steps which are standard in perturbative 
QCD in the momentum subtraction scheme. One usually uses as a subtraction 
point for the three gluon vertex function an Euclidean point with symmetric 
momenta: p\ = p| = p\ = /i 2 . 

In H this non-perturbative minimum subtraction calculation has been per- 
formed, but using an asymmetric subtraction point p\ = p\ = /i 2 ,£>2 = 0. The 
presence of a vanishing momentum induces some subtleties which will be dis- 
cussed later. The running coupling constant computed in || shows a signal of 
perturbative scaling at large scale. 

In this letter we perform the same program with symmetric subtraction points, 
which is the genuine non-perturbative momentum subtraction scheme. We also 
repeat the work done in |J. The whole calculation is achieved with a larger 
statistics, a check of finite volume effects, and a check of scaling when the lattice 
spacing a is varied. 

In section 1 the general principle of the method is recalled and the systematic 
construction of the symmetric momentum points is summarized. In section 2 the 
lattice calculation is described including the checks of color behaviour, of finite 
volume effects, and the discussion of 0(a 2 /i 2 ) and three loop effects. Both the 
scaling in \x and in a are demonstrated. In section 3 we compare our result for 
^JTs w ith other lattice approaches and conclude. 

1 Computing a s from non-perturbative Green 
functions 

In this section we describe the general method used to compute a s and Aq CD in 
the continuum, assuming one is able to compute the Euclidean Green functions of 
QCD in the Landau gauge. The lattice aspect will be treated in the next section. 
The principle of the method is exactly the standard textbook one, generalized 
to non-perturbative QCD ||. 

The Euclidean two point Green function in momentum space writes in the 
Landau Gauge: 

GUT (P, ~P) = G {2) (P 2 )5 aia2 (<W - ^) (1) 
where ai, ai are the color indices ranging from 1 to 8. 
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The three-gluon Green function is equal to the color tensorQ f aia2a3 times a 
momentum dependent function which may be expressed || as a sum of scalar 
functions multiplied by tensors: SA(Pi)P2)Pi)^;MiM2M3' ^ e ah«2M3 being built 
up from ^-|t fc and momenta. In general there is some arbitrariness in choosing 
the tensor basis. One of these may be taken to be the tree level three-gluon 
vertex, projected transversally to the momenta (Landau gauge): 

ir , , , 111, Pi^PifH 



*SS W = IWM - p*U + c y cl - p erm -] II - ( 2 ) 

i=l,3 V Pi I 

while the choice of the other tensors in the tensor basis will be explained below in 
the particular cases. Calling (p\,p\,Pz) the scalar function which multiplies 
the tensor (|2|), the renormalised coupling constant in the considered scheme is 
given by @] 

where 

Z 3 ( M 2 ) = G (2 W (4) 

and fi 2 is the renormalisation scale which will be specified in each scheme. 

The justification of eq. (^|) is standard: the momentum scheme fixes the 
renormalisation constants so that the two-point and three-point renormalised 
Green functions at the renormalisation point take their tree value with the only 
substitution of the bare coupling by the renormalised one. In particular the 
renormalised Gr\p 2 ) takes its tree value, 1/p 2 , at p 2 = fi 2 , which fixes the field 
renormalisation constant (f|) . The renormalised coupling constant is then defined 
so that the three-point Green function is equal to the bare tree level one (with the 
substitution of go by qr) at the symmetric Euclidean point p\ = p\ = p| = fi 2 . 

At the symmetric point only two independent tensors exist in Landau gaugeQ, 
which we choose to be: 



G%« i :r(Pi,P2,P3) = r 



2^ ( Pi ~P2)^(P2 ~P3)/n(P3 

V 2 



(5) 



with T tree defined in (0). To project out G^ 3 ' (/J? , p? , fj?) we contract with the 
appropriate tensor: 

G®{n 2 ,n 2 ,p?)r a2a * = ^GjJ';*" 3 (pi , pz, p 3 ) 



3 In general schemes and gauges the d aia2a3 tensor should also be considered, but not in our 
case as we shall see. 

4 Notice that in the Landau gauge the transversality condition with respect to the external 
momenta, reduces the number of indepcndant tensors as compared to a general covariant gauge. 
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rptree , 



2/i 2 



(6) 



In the following we will call this momentum configuration the "symmetric" one, 
and this defines the MOM scheme. 

We have also considered the MOM scheme defined by subtracting the vertex 
function at the asymmetric Euclidean point p\ = p\ = fi 2 , p 2 = 0. In Landau 
gauge there remains only one tensor, the one in (0), which simplifies: 



Gf 3 A™ 3 (p,0,-p) 



2/ 



110203 



and the scalar factor is extracted via 



^1/^3 



I* 



G (3 V,0,/. 2 



(7) 



(8) 



In the following we will call this momentum configuration the "asymmetric" one. 

Some caution is in order for the latter asymmetric configuration. From ([]]) it 
is clear that 



(p, -p) 



(9) 



for any non vanishing value of the momentum. But when the momentum vanishes, 
the term p^p^/p 2 is undetermined. It could seem quite natural to follow by 
continuity formula (|9|). On the other hand since only the tensor 5^ is defined 
for zero momentum, Gffi(0,0) = 5^ V G^ 2 \Q) leads to replacing the factor 1/3 
by 1/4. Indeed, the Landau gauge condition does not eliminate global gauge 
transformations, and one additional degree of freedom is left at zero momentum. 
This theoretical issue is delicate but it is perfectly obvious that the numerical 
results favor in a dramatic way the factor 1/4. When using the factor 1/3 no sign 
of perturbative scaling can be seen. The factor 1/4 was used in || and we will 
follow the same recipe. 



1.1 Computing A QC d 

The conventional two-loop A* ) is obtained in any scheme from a(/i 2 ) = (?i?(/i 2 ) 2 /(47r) 
by 

where 

da fa 2 /3i 3 fa 4 
a— — = a — t^ol t^ol — ... (11) 

j3 = 11, (3\ = 51 and (3 2 is scheme dependent (P 2 ms = 2857). Integrating 
exactly eq. ( |Tl"P expanded up to order a n+1 , and imposing the asymptotic limit 
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A^ n '/A^(a) — > 1 when a — > 0, leads to the definition 



at two loop and to the three loop A^: 



A(Wl + A^^ A ( C ) (a) (12) 



Pi 



A< 3 > = A< c >(a) 1 + ^+ , 

1 J I 2tt/3o 32tt 2 /? 



x 



exp 



.A \ /VA 
arctan — ; — — arctan — — 



(13) 



2/3 2 VA 

when A = 2/? /3 2 - 4/? 2 > 0. 

One simple criterium has been proposed in to exhibit perturbative scaling: 
when plotting ( [T0|) as a function of \i perturbative scaling implies that A should 
become constant for large enough /i. We will thus try to fit each of the formulae 
(|10|), QI2]) , (pf), expressed in terms of our measured a(/i), as a constant. All these 
formulae converge to the same A when a — > 0. But since our fits are for a in the 
range of 0.3 - 0.5, and since they do not have the same dependence on a they 
should not all fit our data. However, as we shall see, within our errors, acceptable 



fits are possible with fllCf ), and (|13D varying (3 2 on a wide range. This is due to 
the fact that, even with our rather large statistics, the three loop effect, being 
only logarithmic, does not modify strongly enough the variation of A(a) in our 
fitting range although the resulting fitted A depends significantly on the formula 
used and on In other words we have different acceptable fits, with slightly 
different slopes in /i, which lead asymptotically, when a — > 0, to significantly 
different A's. There results a systematic error which cannot be fully eliminated 
until fa is really computed. 

Notwithstanding this problem, we believe that the possibility to fit several of 
these formulae by a constant on a large range of /i and with small statistical errors, 
is an indication that perturbative scaling has been reached. In other words, our 
data show that the uncertainty is of logarithmic type (higher loops), but there is 
no room for significant power corrections. To study power corrections, one has 
to consider with care lower scales [0, and we plan to do that in a forthcoming 
publication [pif] . 



In order to compare different schemes, it is standard to translate the results 
into the MS scheme. Once known in any scheme, a one loop computation is 
enough to yield Aq CD in any other scheme |L0| , [jTTJ] . From JlOf and we get for 
zero flavors 

Amom = 3.334A M g, A ~ M = exp(70/66) ~ 2.888A M g. (14) 
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1.2 Momenta 



In a finite hypercubic volume the momenta are the discrete set of vectors 

y il = lTxn il jL (15) 

where are integer and L is the lattice size. The isometry group for momenta is 
generated by the four reflections p^ — > — p M for \i = 1, • ■ ■ , 4 and the permutations 
between the four directions such as p x <-> p y . Altogether this group has 2 4 4! = 
384 elements. We use fully this symmetry in order to increase the statistics. 
The functions G^ 2 '{p 2 ) and G^'^p 2 , 0,p 2 ) in eqs (ffl) and (|7|) are systematically 
symmetrized over the momenta lying on one given group orbit. The number of 
distinct momenta in an orbit is 384 or a divisor of 384 (when the momentum is 
invariant by some subgroup). 

In the case of G^(p 2 ,p 2 ,p 2 ) in eq. (§J) we furthermore symmetrize over the 
6 permutations of external legs (Bose symmetry). The number of elements in an 
orbit will be 6 x 384 = 2304 or one of its divisors. 

1.3 Triplets of external momenta 

We build all triplets of momenta up to some maximum value of the momentum 
to be specified later. 

In the asymmetric case, there are as many triplets as momenta. For every 
integer number there exists at least one orbit of the isometry group with n^n^ 
equal to that integer. The number of elements is often a much smaller number 
than 384: 1, 8, 16, 24, 64, etc. 

In the symmetric case, n ll n fl has to be an even number: n\ = n 2 , — (n\ +n,2) 2 , 
where the subindices label the external particles, and consequently, n^n^ = 
—n\/2 being an integer, nf is even. It happens that for every even integer, we 
have found at least one orbit. The number of elements in one orbit is often a large 
number, 2304 and 1152 are frequent, 576 and 192 are less. Notwithstanding these 
larger sets, the statistical noise will turn out to be larger in the symmetric case 
than in the asymmetric one. Let us give some examples of symmetric triplets: for 
n 2 = 2: (1100) (0 - 110) (-10 - 10) and its 192 tranformed by the isometry-Bose 
group; for n 2 = 4: (2000) (— 1111) (—1 - 1 - 1 - 1) and its 192 transformed. For 
n 2 = 16: (4000) (-2222) (-2 - 2 - 2 - 2) its 192 transformed. For n 2 = 18 there 
are 6 orbits, for example (4110)(— 1 — 2 — 32)(— 312 — 2) and its 2304 transformed, 
(3300) (-1 - 232) (-2 - 1 - 3 - 2) and its 1152 transformed, etc. 

2 Lattice calculation of a s and Aqcd- 
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2.1 The calculation. 



The calculation has been performed on a QUADRICS QH1, with hypercubic 
lattices of 16 4 and 24 4 sites at /3 = 6.0, 24 4 at /3 = 6.2 ( all three with 1000 
configurations) and 32 4 at f3 — 6.4 (100 configurations), combining the Metropolis 
and the overrelaxation algorithms. 

The configurations have been transformed to the Landau gauge by a com- 
bination of overrelaxation algorithm and Fourier acceleration. We end when 
|<9 M A M | 2 < 10~ 12 and when the spatial integral of A is constant in time to better 
than 10~ 5 . 

We define 

Ujx)-Ul(x) 1 (UJx)-UUx)\ , x 

where fi indicates the unit lattice vector in the direction /x and g is the bare 
coupling constant, and compute the n-point Green functions in momentum space 
from 

• • -Pn) =< Al\{ Pl )A^{ P2 ) ■ ■■A;i( Pn ) > (17) 
where p\ + p 2 + .. + p n = 0, <> indicates the Monte-Carlo average and where 



A;(p) = ^Tr 



M x + A/ 2 ) exp(ip(x + fi/2))\ a 



(18) 



A a being the Gell-Mann matrices and the trace being taken in the 3x3 color 
space. 

We have computed the Fourier transforms up to a maximum momentum of 
~ 3 GeV at (3 = 6, and ~ 4.2 GeV at /3 = 6.2 and f3 = 6.4. These maxima 
correspond to n 2 = n^n M < 16 for (/3, V) = (6.0, 16 4 ) and < 36 for all the other 
cases. 



2.2 Check of the color dependence 

From the color structure of QCD we expect the two point Green functions to be 
proportional to the color tensor 5 aia2 . This is indeed the case to an accuracy of the 
order of 1 %. Furthermore one can prove from gauge symmetry (global and local) 
and Bose symmetry that the three point Green functions have to be proportional 
to f aia ^ in the MOM and MOM schemes. This is indeed the case, but the 
errors now depend on the momentum. For small values of n 2 the agreement is 
of a few percent. The errors increase with n 2 and when n 2 reaches the 30 's the 
error reaches 100 %. This is an indication that the large momenta are grieved 
by noise. Luckily this caution is necessary only for the very few largest values 
of n 2 that we have considered. Indeed we will exclude the points n 2 > 30 from 
our fits. In order to reduce the noise we work from now on with color averaged 
Green functions: ±G^ aia2as f aiCl2Cl3 and |G( 2 ) aia2 <5 aia2 . 
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Finite volume effects ; (3=6.0 
Symmetric case 
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Finite volume effects ; (3=6.0 
Asymmetric case 
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Finite volume effects ; (3=6.0 
Asymmetric case 
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Figure 1: Comparison between the volumes of 16 4 and 24 4 at f3 — 6 for the 
coupling a{fi) (figs, a and b), A$ OM (fig. c) and aJ^L^ (fig. d). No "sinus 
improvement" has been applied here. 
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p 


Volume 


range 


(GeV) 


aA Ms/ ( a V°) s y m 


aA^/(a v /a) asym 


6.0 


16 4 


2.1 


- 3.0 


0.845(25) 


0.809(22) 


6.0 


24 4 


2.1 


- 3.0 


0.811(24) 


0.784(23) 


6.2 


24 4 


2.1 - 


3.85 


0.861(27) 


0.805(21) 


6.4 


32 4 


2.1 - 


3.85 


0.861(45) 


0.793(27) 



Table 1: Ratios aA^/ (a^fa) for different values of (3 after "sinus improvement" 
as explained in section |2.4j . The numbers for a^fa are taken from . The errors 
on aA and a^/a have been combined in quadrature. The scaling invariance is 
particularly striking especially at constant physical volume, i.e. comparing lines 
1,3,4- 



2.3 Finite volume effects 

The finite volume effects can be checked by a comparison of the two calculations 
at /3 = 6.0, Fig [1|. For relatively small //, close to the maximum of a, there is a 
visible decrease of a when the volume is increased. 

For larger /i, the volume dependence is still visible, but reduced to a few 
percent. Comparing the values of A fitted in the asymptotic region given in table 
1, one finds 

A sym (2A)/A sym (lQ) = 0.96 ± 0.02, A asym (2A) / A asym (lQ) = 0.97 ± 0.02 (19) 

which indicates that the finite volume effect affects moderately the asymptotic 
estimate of A. More study is needed to quantify precisely this effect. Still, from 
gross estimates, our largest physical volume, (3 = 6.0, 24 4 , lies presumably within 
5% above the infinite volume limit. 

2.4 Scaling in fi and 0(a 2 p 2 ) effects 

Figs. 0(a,b) show the shape of a(fi). The same shape is seen for the other /3's. 
In fact a(/x) scales in a to a very good accuracy. We keep this study for another 
publication [ |14| 



Turning to the scaling in /i, we see from figs |l|(c,d) that both A's do not 
really show plateaus at large momentum: they go through a maximum around 2 
GeV and fall down later on. Our study shows that this feature cannot be cured 
simply by a three loop effect. Using eq. ( |TB| ) with different values for fi 2 cannot 
lead to acceptable plateaus for all lattice spacings. Since the fall at large /i is 



9 



Results for A with (box) and without (cross) 
^sinus corrections . (3=6.2 , V olume=24 



.40 

.35 - 
.30 

.25 1 

.20 - 

.15 - 

.10 - 



1.0 2.0 3.0 4.0 

Symmetric configuration of momenta 



Results for A with (box) and without (cross) 
sinus corrections . (3=6.2 , Volume=24 




1.0 2.0 3.0 4.0 

Asymmetric configuration of momenta 



Figure 2: The effect of the "sinus improvement" on is illustrated in the case 
of the 1000 configurations at (/?, V) = (6.2, 24 4 ) 
seen in all cases. 



A similar improvement can be 



observed systematically, beyond statistical errors, but decreases when (5 increases, 
we conjecture that we deal with an 0(a 2 p 2 ) effect. 

We have successfully tried a correction which will be described now. It 
starts from the remark that in the lattice Landau gauge, obtained by minimizing 
X^ jX Tr[l — Un(x)\, PpA*(p) does not vanish while p^A^p) does, when A^(p) is 
defined from eq (|18D and where 

a- = 1 < 20 > 

The latter momentum differs from the one in ([T5|) by 0(a 2 p 2 ): p^ ~ p^(l— ^a 2 p^). 
It results that the lattice two point Green function is not really proportional to 
the tensor in ([I]) but to the tensor deduced from (|I|) with p^ substituted by p M , 

We perform a similar change in the tensors used to extract a. The projectors 
in (|D and (§) have been normalized to give 1 when contracted to the tensors which 
multiply G^ 3 ' in (|) and (0) respectively Assuming that the lattice calculations 
is such as to produce the tensors in (|5|) and (|7|) with p^ substituted by p^, there 
would be a bias in our formulae @ and (|j). Indeed the contraction of the 
"tilded" tensors in (|5|) and (|7|) with the tensors in (|6]) and (||) is smaller than one 
and decreases with increasing p. We tentatively correct the bias by dividing the 
result in (|3|) by this factor smaller than one. We shall refer to this as the "sinus 
improvement" . 

For brevity we only show (/?, V) = (6.2, 24 4 ) in fig. |2|. The improvement of 
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p 


6.0 


6.0 


6.2 


6.4 


Volume (lat) 


16 4 


24 4 


24 4 


32 4 


Volume (phys) 


(1.63 fm) 4 


(2.44 fm) 4 


(1.75 fm) 4 


(1.75 fm) 4 


range (GeV) 


2.1 - 3.0 


2.1 - 3.0 


2.1 - 3.85 


2.1 - 3.85 


rVdof 


356(4) 


345(6) 


354(5) 


353(11) 


1.44 


1.12 


1.87 


1.13 


A £ ( MeV ) 

X 2 /dof 


401(5) 


389(7) 


397(6) 


394(14) 


3.13 


1.66 


5.99 


1.43 


A S (MeV) 

X 2 /doi 


314(3) 


303(5) 


313(4) 


312(9) 


1.30 


1.29 


1.02 


1.18 


unimpr. (MeV) 


332(4) 


324(6) 


338(5) 


344(11) 


X 2 /dof 


6.51 


3.71 


12.3 


1.58 



Table 2: Fitted Aj^g for the asymmetric momentum configurations. The "sinus 
improved" A's is used except for the last two lines. To exhibit the a-scaling we 
use the ratios y/a (aA^/ (ay/a)) with o"o = 445 MeV as justified in section pT6| . 
A^ 3 ) has been computed with fa = 1-69/3 2 ms — 4824. The x 2 /dof's correspond 
to the fit in the preceding line. 



the plateaus is dramatic. The large /x fall has been considerably reduced. The 
improvement is confirmed by a reduction of the x 2 P er degree of freedom from 
exceedingly large values to acceptable ones, see tables || and [| 

Of course, this is only an ad hoc 0(a 2 p 2 ) improvement, by no way rigorous 
and systematic. Fitting directly a corrective term of the form 1 — ca 2 p 2 leads 
also to drastically improved \ 2 with best values of c in reasonable agreement 
with the "sinus improvement" (c ~ 1/24). It should be stressed that the sinus 
improvement, and the 1 — ca 2 p 2 fits yield very similar values of A. We may 
thus conclude that the 0(a 2 p 2 ) systematic error on A is moderate after "sinus 
improvement" . 



2.5 Three loop effect. 

A final source of systematic uncertainty comes from our ignorance of fa in the 
MOM scheme. From a preliminary perturbative calculation jT5] in the MOM 
scheme we get fa ~ 4824. On the other hand we unsuccessfully tried to fix fa 
non-perturbatively from our asymptotic fits. The ratio A^/A^, eqs (|10|)- (|T3D , 
drops from 1 when a increases, the drop increasing with fa. As a result, the 
fitted value for A^ 3 ) will decrease as fa increases. Simultaneously the shape of the 
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p 


6.0 


6.0 


6.2 


6.4 


Volume (lat) 


16 4 


24 4 


24 4 


32 4 


Volume (phys) 


(1.63 fm) 4 


(2.44 fm) 4 


(1.75 fm) 4 


(1.75 fm) 4 


range (GeV) 


2.1 - 3.0 


2.1 - 3.0 


2.1 - 3.85 


2.1 - 3.85 


rVdof 


376(6) 


361(6) 


381(8) 


383(20) 


0.38 


1.06 


0.89 


1.50 


A £ ( MeV ) 

X 2 /dof 


442(9) 


425(10) 


444(12) 


455(30) 


3.22 


1.58 


1.89 


1.06 


A S (MeV) 

X 2 /doi 


326(5) 


311(5) 


329(7) 


330(16) 


0.19 


1.25 


1.12 


1.79 


unimpr. (MeV) 


364(6) 


351(7) 


372(8) 


380(21) 


X 2 /dof 


2.99 


1.49 


1.69 


1.24 



Table 3: This table is the analog of table H but for the symmetric momentum 
configurations. 



plateaus are modified. In principle, the requirement of an acceptable x 2 might 
have restricted the admissible domain for f3 2 . Unhappily our preliminary analysis 
did not turn out to be so restrictive. Only for /3 2 below 0.5P 2 ms does the x 2 
become prohibitive, see for example the (3 2 = case (A( 2 )) in tables @ and § It 
might look strange that A^ c ^ fits well while A^ does not, both being two-loop 
formulae. In fact A^ c ^ is only an approximate two- loop formula which can be 
proven to be very close to A^ 3 ) with (3 2 — 8^ ~ 0.66/? 2) ms- 

We therefore cannot do better in the MOM scheme, at present, than to provide 
fits of A^ 3 ) as a function of (3 2 . For comparison we also provide the same analysis 
in the MOM scheme. The maximum value for f3 2 which we consider is f3 2 = 
2/^2 Ms = 5714 since, for such a large value, the term oc a 4 in the (3 function 
(|Tl"D is of the same order as the term oc a 3 for our range of a. If (3 2 was larger 
than that, the perturbative expansion would be dubious, and the evidence for 
perturbative scaling shown by our data would appear as a miraculous fake. 

2.6 Scaling in a 

The "sinus improved" A's exhibit a very clear scaling when /x varies above 2.1 
GeV, as can be seen from the quality of the plateaus in figs || and |] and from the 
X 2 per d.o.f. In this subsection we want to study further the scaling when (3, i.e. 
a, is varied. 

Since A depends linearly on a -1 , the consistency of our fits can only be checked 
through a spacing independent ratio. We use the ratios a A/ '(ay/a), see table [3], 
where ay/a is the central value of string tension computed in ||. 
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In order to write A in physical units we then multiply all ratios by one global 
scale factor: yfa = 445 MeV tuned to the central value of a very recent fit [13 



from the K* mass: a (/3 = 6.2) = 2.75(18) GeV. We take the central value: 
a' 1 (/3 = 6.2) = 2.75 GeV, whence a" 1 ^ = 6.0) = 1.966 GeV and a" 1 ^ = 6.4) = 
3.664 GeV. 

This leads to the plots in fig. |3|. The presence of nice plateaus is striking. We 
fit the average A on these plateaus, for scales never smaller than 2.1 GeV, and 
as high as allowed by lattice effects. The results are presented in tables [I], |2| and 
H The fits for A^ c ^ and for a large range of A^ yield a x 2 P er degree of freedom 
smaller than 1.5. Scaling in the lattice spacing is striking, especially for those 
lattice parameters which correspond to a similar physical volume of ~ (1.7 fm) 4 , 
i.e. (/3,V) = (6.0, 16 4 ), (6.2, 24 4 ) and (6.4, 32 4 ). They average to: 



A^L = 378(6) MeV (symmetric) A^ = 355(4) MeV (asymmetric) 



A S( , 2 = 1MPaa = 4824) = { W (s = ic) (21) 

where the errors are only statistical. {3% = 4824 results from our preliminary 
calculation [[DJ in the asymmetric scheme. For comparison we provide the result 
with the same fa i n the symmetric case. The result at the larger volume of 
~ (2.44 fm) 4 , (/3,V) = (6.0, 24 4 ), presumably close to the infinite volume limit 
(section |2.3|) , is: 

A^L = 361(6) MeV (symmetric) A^l = 345(6) MeV (asymmetric) 

a (3) to _ i fiq a _ _ 4 o 94 n _ / 311(5) MeV (symmetric) . , 

a ms ^ 2 - i ' Dy ^2MS - ~ \ 3 3( 5 ) MeV ( asymme t r i c ) W 

Varying (3 2 we find acceptable x 2 ' s from /9 2 ~ 0.5/9 2 j^g = 1428 up to beyond 
2/^2 Ms = 5714 which we take as the maximum perturbatively consistent value, 
see section 12.51. In this range of fo the fitted A^ 3 ) have, to a surprisingly good 



approximation, a linear dependence on We provide the result in the next 
section. 

Finally it is worth mentioning that we have also checked scaling of a in a 
over the whole range in fi, including the small values. We leave this point for a 
forthcoming publication |14 



3 Discussions and conclusions 

There is scaling, as can be seen first from the plateaus of A as a function of 
the momentum scale //, and second from the striking agreement of the runs for 
different /?'s. We now quote our final results from our largest physical volume, 



13 



Symmetric case ; (3=6.0 ; Volume=16 




Asymmetric case ; (3=6.0 ; Volume=16 
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Symmetric case ; (3=6.4 ; Volume=32 
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Figure 3: The fits for aA^/(a^)^a (with ^ = 445 MeV), including the 



(fa V) = (6.0, 24 4 ), which we estimate to give values of A less than 5% from the 
infinite volume limit. 

The analysis for symmetric momentum configurations is better grounded the- 
oretically since it avoids the delicate problem of zero momentum. On the other 
hand, this analysis is noisier than the asymmetric one which exhibits beautiful 
plateaus. The good agreement of these two analyses allows a sort of reciprocal 
support. 

Several other lattice estimates of A have been performed. The ALPHA col- 
laboration, 0, quotes Aj^g = 251(21) MeV. Other results are 244(8) MeV §, 
293(18)| MeV §], 340(50) §. 

Our results for A happen to be very sensitive to the three loop effect but fa 
cannot be fitted non perturbatively from our data. A wide range, 0.5/3 2 y^ = 
1428 < fa < 2/3 2 j^g = 5714 is allowed, in which the three loop A^ 3 ^ can be 
approximated by the following formulae: 



ASL = f 412 - 59-^- ± (» ! I\ [<V 

2, MS 



MS \ "~ ""ft 

^2. MS 



1 (/? = 6.0) , 

(symmetric) 



A^L = I 382 - 46-^- ± 5 I MeV 

^MS / 



MS \— 



1.97GeV 

a~\p = 6.0) 
1.97GeV 



(asymmetric) (23) 



Comparing the results in both schemes seems to indicate that the fa's in 
MOM and MOM schemes are not too different. A calculation of fa in the MOM 
scheme would be most welcome. 

Our preliminary computation of fa in the MOM scheme, [JI5] , uses the results 
of |16| and yields a value of /? 2 mom ~ ^ Ms — 4824. Our final result is then 



A^=(303±5MeV) a i 1( 9 /3 7Ge ^° ) (asymmetric) (24) 
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